
#ifndef __GSL_EXT_MATRIX_H__

#include <gsl/gsl_matrix.h>

double 
gsl_ext_matrix_get_max_abs_diag_element(const gsl_matrix *A);

double 
gsl_ext_matrix_get_max_abs_offdiag_element(const gsl_matrix *A);

void 
gsl_ext_matrix_compute_gerschgorin_bounds(const gsl_matrix *A, 
                                          double *minev, 
                                          double *maxev);

#define __GSL_EXT_MATRIX_H__

#endif
